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ABSTRACT 

Calculations are performed to demonstrate the deviations from Gaussian that occur in the spectral line 
profiles of a linear maser as a result of the amplification process. Near-Gaussian profiles are presented 
for bright, interstellar 22 GHz water masers obtained from high resolution Very Long Baseline Array 
(VLBA) observations of W3 IRS 5. For the profiles to be so close to Gaussian, the calculations indicate 
that these masers must originate in quite hot gas with temperatures greater than 1200 K — a conclusion 
that is supportive of C-type shocks as the origin of these masers. In addition, the degree of saturation of 
these masers must be less than approximately one-third, from which it follows that the beaming angles 
are less than about 10 ster and the actual luminosities are modest. If spectral profiles that are as close 
to Gaussian as the profiles presented in this initial investigation are found to occur widely, they can be 
valuable diagnostics for the environments of astrophysical masers. 

Subject headings: shock waves — line: profiles 
radio lines: ISM 
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1. INTRODUCTION 

With increasing intensity, the spectral line profiles that 
are calculated for astrophysical masers first become nar- 
rower and then rebroaden when the masers become satu- 
rated. At high saturation, the profile in Doppler velocity 
becomes the same as the Maxwellian (Gaussian) velocity 
distribution of the masing molecules. This behavior has 
been established to be independent of whether the masers 
are treated as linear (Litvak 1971), or as having finite cross 
sections such as disks and spheres (Emmering & Watson 
1994; Wyld & Watson 2002, in preparation). It seems to 
be less widely recognized that the profiles also can be negli- 
gibly different from Gaussians at low intensities where the 
maser is unsaturated. We demonstrate how the deviations 
of the profile from Gaussian can (along with the breadth of 
the profile) then be used to obtain information about the 
masers. In particular, lack of knowledge about the solid 
angle into which the maser radiation is beamed — which is 
necessary to obtain the luminosity and the energy density 
from the observed radiative flux — has been an outstand- 
ing weakness in the interpretation of astrophysical masers. 

At low angular resolution, the observed spectral line pro- 
files are not expected to be indicative of the profiles of in- 
dividual masers because multiple masing components may 
be present within a telescope beam. However, the pro- 
files presented here have been observed at very high angu- 
lar resolution with the Very Long Baseline Array (VLBA) 
and are quite close to Gaussian. Most likely, the devia- 
tions from Gaussian in these profiles are upper bounds to 
the actual deviations in the profiles of individual masers. 

2. OBSERVATIONS OF NEAR-GAUSSIAN MASER PROFILES 

The 22 GHz H 2 maser profiles toward W3 IRS 5 shown 
in Figure 1 were obtained with the VLBA in order to mea- 
sure the circular polarization of the maser radiation, and 
hence to infer the strengths of the magnetic fields. Also 
shown in this Figure are the Gaussian fits and the differ- 
ences between the fits and the data. The observational 



procedures to obtain these profiles have been described 
in a previous Letter (Sarma, Troland, & Romney 2001). 
Here we examine the near-Gaussian nature of the profiles. 
W3 IRS 5 is an infrared double source located between 
the H II regions A and B in the W3 giant molecular cloud 
(e.g., Tieftrunk et al. 1997). Claussen et al. (1994) found 
seven hypercompact continuum sources (diameter < 700 
AU: Tieftrunk et al. 1997) in the W3 IRS 5 region with the 
Very Large Array (VLA). They interpreted these sources 
as stellar winds from embedded B stars (cf. Tieftrunk et al. 
1997). These authors also found that some of the brightest 
H2O masers coincide with the locations of 15 GHz contin- 
uum images of the hypercompact sources. From the flux 
in the observed profiles and the maser sizes (assumed un- 
resolved; hence sizes are the half-power beamwidth, which 
is ^0.7 mas), we find that the minimum brightness tem- 
peratures Tb for the masers in Figure 1 are 1 to 4 x 10 12 
K. 

3. CALCULATIONS 

The equation of transfer for the dimensionless intensity 
I(v) of a ray of maser radiation at Doppler velocity v can 
be expressed as a function of the optical depth parameter 
s in the linear maser idealization as 



dl 
ds 



I 



l + (l + g u /gi)I 



exp (-v 2 /0.36v 2 ) (1) 



Here, the I = I(v) is in units of the saturation in- 
tensity I s — (8irhis 3 r /c 3 g u AAfl), where v is the fre- 
quency of the transition, T is the "phenomenological" de- 
cay rate of the molecular state, g u ,9i are the degenera- 
cies of the upper and lower states, A is the Einstein A- 
valuc, and AQ is the beaming angle for the maser ra- 
diation. Also, the source function due to spontaneous 
emission S = {(g u AAQ/ 'AnT)(A/ 'AA)] where A is the phe- 
nomenological pumping rate into the upper state and AA 
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Fig. 1. — VLB A 22 GHz H2O maser profiles toward W3 IRS 5. The crosses represent the observed data, and the smooth curves represent 
the fitted Gaussians. The dots in the lower panel represent the residuals from the Gaussian fit that have been estimated at each observed 
velocity. The error bar along the intensity axis at each observed point (as measured from the rms in the images of Stokes V = right minus 
left circular polarization) is 16 mjy beam -1 . The lincwidths (FWHM) are also given in the figure. Other details of the observations are given 
in Sarma ct al. (2001). The maser labels used are also taken from that paper. Note that the residuals for the masers NW and SE have been 
multiplied by factors of 5 and 4 respectively to increase their visibility. 



is the difference between the pumping rates into the upper 
and lower states. The dispersion (FWHM) in the thermal 
velocities of the masing molecules is v t . If spontaneous 
emission is ignored, equation 1 can be integrated to yield 
Iexp[(l+g u /gi)I} = I c exp[sexp(-u 2 /0.36uf)] for I c < I, 
where I c is the intensity of the continuum radiation that 
is incident at the far end of the maser which serves as the 
"seed" radiation for the maser. The I(v) that is found by 
integrating equation 1 is fit to a Gaussian using the least 
squares procedure. The result 5 of dividing the integral of 
the squares of the normalized deviations by the linewidth 
Av 1/2 (FWHM) of the computed profile 

5 = {J [I(v) - ai eM-v 2 /a 2 )} 2 dv}/I%Av 1/2 (2) 

is then adopted as a quantitative measure of the deviation 
of the spectral line profile from a Gaussian. Here, a\ and 
a 2 are the parameters of the Gaussian that are obtained 
in the least squares fitting and I p = I(v = 0) is the peak 
intensity within the spectral line. Note that <5 is indepen- 
dent of v t (that is, of the kinetic temperature of the gas) 
and that Avi/ 2 — v t x f(I c ,I p ) when S is ignored. The 
function f(I c ,I p ) — ► 1 as I p (which also is a measure of 
the degree of saturation) becomes much greater than one. 
In Figure 2 examples of the variations of the calculated 5 
and Avi/2 are shown as a function of the peak intensity I p 
(or the degree of saturation) for I c = 10~ 9 , which we will 
reason subsequently is representative for the data in § 2. 
Clearly, the main deviations in I(v) from a Gaussian occur 
in the region of partial saturation. The two curves for the 
linewidth have different values for vt, and are thus shifted 
vertically from one another by a constant factor. The hori- 
zontal line in Figure 2 intersects these two curves at differ- 
ent points, and is used to indicate that a specific linewidth 
ordinarily is compatible with a range of gas temperatures 



and a range of degrees of saturation. The intersections 
determine the corresponding I p and 5 when a linewidth of 
0.5 results from gas at the two temperatures. In Figure 3, 
the relationship is shown between the intensity I p and the 
gas temperature (scaled) which yield a specific linewidth 
Avi/2, for several values of I c . Note that the information 
in Figure 3 can be scaled as described in the caption for 
any linewidth and molecular species. The "kurtosis" — a 
quantitative measure of whether the profile is more [pos- 
itive] or less [negative] peaked than the Gaussian at line 
center (e.g., Press, Flannery, & Teukolsky 1986) — also 
is shown in Figure 2. That the kurtosis changes sign at 
essentially the value of I p where the linewidth is a mini- 
mum and where the saturation begins to be significant is 
a further indication of the change in the character of the 
spectral line at this point. 

We find that, at least for the ranges that are consid- 
ered here for the parameters, spontaneous emission yields 
results that differ negligibly from those obtained for inci- 
dent continuum radiation with I c = S. We thus present 
results labelled only according to the I c , but with the un- 
derstanding that they are applicable for S — I c as well. 

4. INTERPRETATION OF THE 22 GHZ H 2 MASER 
PROFILES 

The relevant characteristics for the profiles in Figure 1 
are (Av 1/2 km s"\ 5) = C (0.7, 2.4 x 10~ 3 ), NW (0.7, 
2.4 x 10~ 4 ), SE (0.5, 1.8 x 10~ 4 ), and for a fourth profile 
observed by Sarma et al. (2001) but not shown in the fig- 
ure, (0.5, 8.5 x 10~ 4 ). At least within the relevant range 
of values for I c that is considered in Figure 3, the values 
of I p at which specific S occur are essentially independent 
of I c when 5 < 10~ 3 . Hence, to be compatible with the 
observed profiles in § 2 for which three have S < 10~ 3 , the 
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Fig. 2. — The lower panel shows the representative lincwidths 
At>i/2 (arbitrary units) that are calculated as a function of the log 
of the peak intensity I p for masing gas at two gas temperatures. The 
horizontal line at Aui/2 = 0.5 represents an observed linewidth. Its 
intersections with the curves for the calculated linewidths determine 
the gas temperatures and intensities that arc compatible with this 
linewidth. The deviation of the profile from a Gaussian as measured 
by <5 is given by the dotted curve on which the crosses mark 8 = 3 
X 10 ~ 4 that is representative for the observed profiles in Figure 1. 
Though the results of calculations arc presented in this figure only 
for Ic = 10~ 9 , they are indicative for the range of choices for I p that 
are relevant here (i.e., those in Figure 3). The upper panel shows 
the kurtosis for these profiles. Note that 8 and the kurtosis do not 
depend on the temperature of the gas. 



intensity must either be quite large (> 100) or it must be 
somewhat less than 1 according to Figure 2. We reject the 
high saturation solutions because (i) the required satura- 
tion seems excessive, but more objectively because (ii) the 
gas temperatures T that would be indicated (~ 100 K for 
Aui/2 = 0.5 km s~ 4 ) in Figure 3 are below what seems 
to be required for the chemistry and pumping of these 
masers, and (iii) the 22 GHz H2O maser feature probably 
consists of a blend of three hyperfine components. Hyper- 
fine structure is omitted here, but would cause irregular 
profiles at temperatures of a few hundred K. At these tem- 
peratures, v t is too small in comparison with the hyperfine 
splittings for the components to merge completely. In ad- 
dition, (iv) the "kurtosis" of the four line profiles that are 
observed is positive as is expected from Figure 2 for maser 
profiles at I p < 0.3, but not for larger I p . Also, the resid- 
uals at the intensity peak for two of the three profiles in 
Figure 1 (and for a fourth profile that is not shown) are 
positive, directly indicating that the observed profiles are 
more sharply peaked (positive kurtosis) than the Gaussian 
fit. For the third profile shown in Figure 1 (maser SE), the 
peak in the residual is shifted by one velocity channel from 
the peak in the observed profile. 

Now consider the likely ranges for I c ~ 10~ 9 T c Af2 and 
£ ~ 10~ 9 Afi(A/AA), where the background continuum 
intensity has been expressed as a brightness temperature 
T c and the usual estimate for the 22 GHz masers T = 2 
s~ 4 has been adopted. In calculations for the pumping of 
these masers, the fractional inversion AA/A typically is a 
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Fig. 3. — The scaled gas temperature and the log of the peak 
intensity I p which yield a spectral linewidth Av 1 / 2 (FWHM) in 
calculations analogous to those for Figure 2. The scaled gas tem- 
peratures must be multiplied by the factor t =(Ai>i/2/0.5 km 
s -1 ) 2 X (M/Mh 2 o) to convert them to actual gas temperatures in 
Kelvins, where M is the mass of a molecule of the masing species. 
The various curves (solid lines) arc computed for different choices of 
I c = lfj- 7 , 10- 8 , 10- 9 , 10- 1( », 10- 11 as indicated. The horizontal 
dashed line at log I p = —0.5 marks the approximate upper limit in 
Figure 2 for the near-Gaussian profiles for which <5 < 3 X 10~ 4 . The 
dotted line represents I c = 10~ S I P , which follows from the lower 
limit to the observed brightness temperatures for the masers in § 
2. According to the interpretation, the allowed gas temperatures 
must be on the high temperature side of this line as indicated by the 
arrows. 



few percent (e.g., Anderson & Watson 1993) to give S ~ 
10 _7 A£L Plausible values for T c range from little more 
than the 3K cosmic background to the 10 4 K that is char- 
acteristic of H II regions. Then, since it is unimportant 
for the calculated spectra whether the seed radiation is 
provided by spontaneous or by background radiation, for 
convenience we can express the likely range for the seed 
radiation in terms of I c alone with 10 4 > T c (K) > 10 2 . 
Because I p can be related to 5, it is useful to eliminate AQ 
and express I c as I p T c /T B = 10-%[(T C /10 4 K)/(T B /10 12 
K)], or I c < 10~ 8 ip here, since T B = 10 12 K is a lower limit 
for the data in § 2. The locations on the curves in Figure 
3 where I c — 10~ 8 I p are indicated by the dotted line. The 
upper end of this line is at I p = 0.3 (and I c = 3 x 10 ) — 
the value at which S = 3 x 10~ 4 in Figure 1. Hence, the 
allowed region in Figure 3 is that indicated by the arrows, 
and this dotted line indicates the lowest gas temperatures 
that are compatible with the two profiles in § 2 that are 
closest to Gaussian profiles — and probably with the other 
two profiles as well in view of the sense of the uncertainties 
in T c , Tb and 5. In turn, AQ can be related to the inten- 
sities and brightness temperatures of the masers through 
An = (4whuT/ckg u A)I p /T B ~ 3 x 10-%/(T B /10 12 K) 
ster. 

It is evident that the foregoing analysis indicates that 
the maser profiles described in § 2 originate in gas that is 
quite hot — if taken at face value, in gas with tempera- 
tures greater than about 1200 K and 2500 K based on the 
spectral lines with widths of 0.5 (S = 1.8 x 10~ 4 , and 8.5 
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x 1CT 4 ) and 0.7 km s" 1 (S = 2.4 x 1CT 4 ), respectively. 

For simplicity, the foregoing analysis has been based on 
a two-level maser transition and the linear maser ideal- 
ization. The 22 GHz transition probably is a blend of 
three hypcrfinc components. Though this can influence 
the analysis in detail, our initial calculations indicate that 
the basic conclusion that the observations imply I p < 0.3 
and high gas temperatures greater than about 1000 K will 
be unaltered. Perhaps surprisingly, the line profiles that 
result from the three merged hypcrfine components still 
tend to be good Gaussians except at the lowest temper- 
atures. Again for simplicity, we have not dealt explicitly 
with the possibility that the rate r„ for the relaxation of 
molecular velocities by the trapping of infrared radiation 
can be greater than the decay rate T of the state (Goldrcich 
& Kwan 1974). If so, will replace T in our calculations 
(Ncdoluha & Watson 1991). However, T will be reduced 
by the trapping, and r„ will have approximately the rate 
(2 s _1 ) that we have used in the foregoing discussion (An- 
derson & Watson 1993). Any changes in our conclusions 
are expected to be minimal. The linear maser idealization 
is considered to be applicable because astrophysical masers 
(especially the 22 GHz masers) are thought to be highly 
elongated with very small beaming angles so that the rays 
are nearly parallel. The A£l < 10~ 4 ster inferred above 
supports this premise for the masers of § 2. We cannot, of 
course, exclude the possibility that velocity and/or excita- 
tion gradients are present within the maser and have just 
the right variation to enhance the lincwidth while main- 
taining a near-Gaussian shape for the overall profile. It 
seems unlikely that a velocity gradient which is simply 
constant along the line of sight will do this. A velocity 
gradient will tend to flatten the peak and broaden the 
overall profile. When I p > 0.3 in Figure 2, the peak of the 
profile already is more flat than a Gaussian. A constant 
velocity gradient in this regime would only increase the 
deviation from a Gaussian. When J p <1 and the maser 
is unsaturated, a velocity gradient that is large enough to 
increase the linewidth significantly would be expected to 
cause excessive distortion in the line profle. The remaining 
possibility in a constant velocity gradient is that velocity 
relaxation is significant and that T v is greater than the 
rate for stimulated emission, which in turn is greater than 
T. This would also correspond to the regime I p < 0.3 



in Figure 2. Although the peak of the profile does not 
seem to be flattened in this case, the available calculations 
(Ncdoluha & Watson 1988) also suggest that the velocity 
gradient causes little increase in the linewidth — at least 
near the minimum for the linewidth in Figure 2. Our in- 
terpretation of the data in Figure 1 would then be largely 
unchanged. Further, detailed calculations as a function of 
the relevant parameters are desirable, however. 

5. DISCUSSION 

An analysis of the line strengths of the masing transi- 
tions of water at higher frequencies by Melnick et al. (1993) 
has indicated that these masers originate in a gas that is 
quite hot (> 900 K). Although 22 GHz H 2 masers do 
occur along the same lines of sight, the possibility that 
they are at different locations along those lines and might 
originate in gas at the lower temperatures could not be ex- 
cluded. Our interpretation of the Gaussian nature of the 
22 GHz profiles indicates that strong 22 GHz H2O masers 
do originate in gas at similarly high temperatures. Shock 
waves have long been recognized as likely environments for 
interstellar water masers (Shmeld, Strclnitskii, & Muzylev 
1976). Recent attention has been focussed on the specific 
properties of J-type shocks and C-type shocks. In the for- 
mer, the 22 GHz masing has been found to originate in 
a relatively narrow range of temperatures near 400 K at 
which hydrogen molecules recombine. In contrast, 22 GHz 
masing at gas temperatures up to 3000 K and higher can 
occur in C-type shocks (Kaufman & Neufeld 1996). Ev- 
idently, our analysis is compatible with an origin for the 
masers in C-type shocks, but not J-type shocks. Another 
implication of the relatively low degree of saturation of 
these masers is that the standard Zceman interpretation of 
their observed circular polarization (Watson & Wyld 2001) 
is more likely to be valid (Sarma et al. 2001). Clearly, more 
extensive observational efforts will be necessary to deter- 
mine whether 22 GHz masers, and other masing species, 
commonly have near-Gaussian spectral profiles which can 
be used as diagnostic tools. 
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